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We present a self-consistent approach to describe ambipolar tunneling in asymmetrical double 
quantum wells under steady-state excitation and extend the results to the case of tunneling from a 
near-surface quantum well to surface states. The results of the model compare very well with the 
behavior observed in photoluminescence experiments in InGaAs/InP asymmetric double quantum 
wells and in near-surface AlGaAs /GaAs single quantum wells. 
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I. INTRODUCTION 



Tunneling of electrons and holes in quasi-two-dimensional semiconductor heterostructures is the object of large 
interest both for its physical interest - it is one of the most important quantum-mechanical effect observed in low- 
dimensional structures - and for its role in several nanometric devices and applications. Much work has been done 
concerning tunneling in symmetrical and asymmetrical double quantum wells ( ADQW) . 

Tunneling from a quantum well (QW) to surface states has received instead little attention, despite the great 
importance of gaining control over this mechanism. For a quantum well built in the neighborhood of an unpassivated 
surface, an extra non-radiative recombination channel becomes available to electrons and holes if they can tunnel to 
surface states, with consequent loss in emission efficiency from the quantum well. The importance of this effect has 
been demonstrated experimentally in recent papers and its dependence on surface- barrier thickness investigated [Q-Q . 

In the interpretation of the experimental results, both tunneling to surface states and tunneling between ADQW, 
a problem arises. Although in some cases hole-tunneling rates can be comparable to the electronic ones, e.g. when 
heavy holes in one well move to a light hole state in the other well |^ , direct-gap III-V semiconductors have in general 
rather different tunneling probabilities for electrons and holes due to the different effective masses . The different 
tunneling probability for the two carriers causes a dipole electric field across the tunneling barrier to develop [p|-pd|, 
so as to induce - via the quantum confined Stark effect - a peak shift of the excitonic recombination, and to affect 
the tunneling probabilities towards an ambipolar regime, with equal tunneling currents for electrons and holes. 

The appearance of an electric field due to the spatial separation of electrons and holes in ADQW systems was 
already discussed theoretically p3| in the framework of exciton tunneling under impulsive excitation. In this paper we 
study the tunneling of unpaired electrons and holes both in ADQW and in near-surface QW systems under steady- 
state excitation. Unlike the case of impulsive excitation, not discussed in this paper, the ambipolar regime is shown to 
be reached in this case for any excitation intensity. The consequently modified tunneling properties are successfully 
compared with available experiments, i.e. dependence of the emission efficiency ratio in ADQW systems on the 
excitation power and the tunneling-barrier width j^] , dependence of the tunneling current on the barrier thickness , 
and dependence of the Stark shift on the excitation level in a near-surface QW |q| . 



II. TUNNELING BETWEEN TWO ASYMMETRIC QUANTUM WELLS 

Here we describe the steady-state photoluminescence from two coupled asymmetric quantum wells under constant 
irradiation. The charges photogenerated inside the wells relax almost instantaneously with respect to the other 
relevant time scales to the lowest band of the respective well and fill it according to the Pauli principle. Electron-hole 
interaction leads to exciton formation. Tunneling between the two wells is essentially restricted to unpaired electrons 
and holes when the tunneling process is a small perturbation for the nearly uncoupled wells: as a matter of fact, the 
tunneling current densities correspond to a first-order process for unpaired electrons and holes and to a second-order 
process for excitons. On the other hand photoluminescence is restricted only to excitonic recombination. 

Let oi and 02 be the widths of the two quantum wells and b the width of the barrier between them. We suppose 

(2) (2) 

fli > 02 so that the bottom of the el and hhl bands of well 2, E^^ and Ef^j^^, are higher in energy than those of well 
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1, E^]^ and -Ej^^j^- Let Gi and G2 be the generation current densities of electron-hole pairs in the two wells. If Ui and 
Pi are the steady-state concentrations (number of particles per unit area) of electrons and holes in the well i = 1,2, 
the following rate equations hold: 

— Gi + Je — Xinipi (la) 
= Gi + Jh- Xifiipi (lb) 

= G2 - Je - \2n2P2 (Ic) 

= G2 - J/i - X2n2P2 (Id) 
= AiTiipi - Ii (le) 

= A2 712^2-^2- (If) 

The first two couples of equations are the rate equations for the concentrations of unpaired electrons and holes 
in well 1 and 2, respectively, and contain, in the order, the photogeneration and tunneling current densities of the 
unpaired charges and the generation current density of excitons. The last two equations are the rate equations for 
the exciton concentrations in well 1 and 2 and contain the photoluminescence current densities Ii and I2 which are 
the quantities we want to evaluate and compare with experimental results. 

The generation current density of excitons in well i is assumed proportional to the concentrations of the unpaired 
charges in the same well, XiUiPi, where Xi is the bimolecular exciton formation coefhcient possibly dependent on the 
well size Q. 

A little more complicated analysis is needed for the tunneling current densities Je and J/j. Transfer of electrons 
(holes) from the narrower well to the larger one is realized in a non-coherent two-step process. Quantum coherent 
tunneling of electrons (holes) from an occupied state of the el (hhl) band of well 2 to an equal-energy empty state 
of the el (hhl) band of well 1 is followed by thermalization toward the lower energy states of well 1. When the 
barrier width b is not too small, which is also the range of validity of the above rate equations, the quantum coherent 
tunneling process characterized by a time growing exponentially with b gets much slower than the thermalization 
process in well 1 (^ 1 ps [^) and this last can be neglected. 

The current densities for quantum coherent tunneling are approximately proportional to the charge concentrations 
in well 2 and the proportionality factor, namely the tunneling rate, is quite different for electrons and holes due to 
their different effective masses. Therefore, in a steady-state situation when J,, — Jh, the concentrations of electrons 
and holes in each single well must be different and an electric field gets established in the barrier between the two 
wells. Moreover the electric field modifies the electron and hole tunneling rates. The direction of the field is simply 
understood from the values of the tunneling rates at zero field. As electron tunneling rate is expected in this case to 
be much greater than the hole tunneling rate, electrons accumulate at well 1 and in a steady-state situation ni > pi- 
The electric field is directed from well 2 to well 1 and its value is given by 

£0£r 

where n = rii — pi = P2 — n2 and Sr is the permittivity of the barrier material. 

The tunneling current densities depend on the electric field F through the dependence of the charge concentrations 
on F as well as through the dependence of the tunneling rate on the energy shift between the bands of the two wells 
induced by F. Assuming that the electric field is completely shielded inside the wells, this energy shift amounts to 
eFb. For a given value of F , the tunneling current densities can be evaluated within perturbation theory in the case 
of a barrier not too thin. At the first order the transition rate induced by a constant perturbation V between initial 
and final states continuously distributed in energy with densities dNi/de and dNf/de, respectively, is [ p^ 

^\{^,e\V\f,e')f 6ie- e')dNf{e')dN,ie) ^ J ^\{^,e\V\f,e)\' ^{e)^ie) de. (3) 

This formula, when applied to the electron and hole states of the two uncoupled wells and divided by the transverse 
area A of the heterostructure, gives the following expressions for the tunneling current densities 
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The integration takes into account the transitions from aU the occupied states in the el {hhl) band of weU 2 (initial 
states) to the empty states in the el (hhl) band of weU 1 (final states) at the same energy e measured from the 
bottom of el {hhl) band of well 2. The densities of the available initial (final) states are obtained by multiplying the 
number of states per unit energy A v^'^^ [A v^^^) with the appropriate occupation probability. Here i/e*'' = mi^ j-Kfi} 
f l-Kh^), i=l,2, are the densities per unit area and unit energy of a two dimensional ideal gas of fermions 
with mass rvt^e' (jn^^^ and j{x) = l/[exp(a;) + 1] is the Fermi function. The Fermi energies for electrons (holes) are 
related to the concentrations in the corresponding well: e^-* — nh^ rii / Tm!i^ +-^ei'' +&Fb — and e^-* = -kU^ n2 / m'e^ 

T^fi^Vil^h')- Finally, the perturbation potential V^, {Vh) is 
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the potential of the heterostructure (& finite) coupling the electron (hole) states evaluated by considering the wells 
uncoupled (& infinite). Due to the exponentially vanishing tails of the electron (hole) wavef unctions, the relevant 
contribution to the matrix elements comes only from the barrier region between the two wells p^ . In this region 
the potential Ve (T4) has magnitude of the order of the conduction (valence) band offset AE'c (AEv). More accurate 
expressions for Ve (Vh), e.g. taking into account the distortion due to the electric field, have little influence on the 
final result and will be neglected. Explicit expressions of the electron and hole states $ and analytical evaluation of 
the matrix elements are given in Appendix A. 

The formulas for the tunneling current densities can be simplified by noting that, due to the exponentially large 
difference between electron and hole tunneling rates, we expect n2 <C ni. As a consequence, the electron Fermi energy 
i.e. the effective region of integration in Eq. (^), is very small. By approximating all factors to be integrated in 
(0) but the Fermi function in well 2 with their value at e = (the bottom of the el band of well 2) we get Je = n2/Te, 

where n2 — J v^^^ ] ^(e — e^p)/kBT^ de and 
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In the case of holes we expect p2 ^ Pi so that e^'' is not small. However, in this case < (now the zero of energy 
is the bottom of the hhl band of well 2) and /((e — e^P)/kBT) ~ 0. Neglecting the smooth dependence of the matrix 
element on e in (H) we get Jh — Vil'^h where P2 = J f ((^ ~ ^^f^) /ksT^ de and 
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Note that both the tunneling rates I/te and l/r/j depend on the electric field through the matrix elements (see 
Appendix A). In addition the electron tunneling rate depends on F through the Fermi energy of the well 1. 

If we think that the electric field, i.e. n, and the tunneling rates are known, we get the following solution for the 
initial set of rate equations 
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This result allows us to find the steady-state values of the electric field, of the charge concentrations and of the 
tunneling rates by a recursive method. We proceed in the following manner. First of all we fix an electric field value 
corresponding to some negative electric charge concentration n in the well 1. Then we evaluate the hole tunneling 
rate and give a starting value to the electron tunneling rate (for instance that obtained for ni = n) so we can deduce 
some charge concentrations by Eq. (^). From these concentrations we evaluate a new electron tunneling rate by 
Eq. and from that new charge concentrations. After few iterations, the electron tunneling rate and the charge 
concentrations converge to the solution corresponding to the value of the electric field fixed at the beginning. At this 
point we compare the values of the tunneling current densities. If Je > Jh the electric field is increased, if Je < Jh 
the electric field is decreased. The procedure is repeated until the steady-state condition Jg = Jh is reached. Note 
that the existence of the ambipolar regime does not rely upon the value of the excitation intensity, i.e. Gi and G2. 
However, the value of the electric field established in the steady state strongly depends on the excitation intensity. 

We illustrate the behavior of the model by considering a practical example close to the experimental situation 
investigated by Sauer, Thonke and Tsang We consider two /rto.53Gao.47^s quantum wells imbedded between InP 
barriers with ai = 100 A and 02 = 60 A. We use the following material parameters at T = 4.2 K AEc = 0.195 

eV, AEy = 0.293 eV, mi^^ = mi^' = 0.044 to, to^^^' = mj^^^ = 0.38 to, to being the free electron mass, and Er — 13.9. 
Moreover we put Ai = A2 = 6 cm'^ s~^ |lj]. Since aai, aa2 ^ 1, where a is the optical absorption coefficient for the 
pump light, we take G1/G2 — 0,1/0,2 and G2 = P/hv where P is the absorbed power density and hv — 2.41 eV is the 
photon energy corresponding to a 514-nm laser light. Note that the transverse area A shown in Eq. (||) and Eq. (|^) 
cancels out with the A"^ from the squared matrix elements and is an irrelevant parameter. 

In Fig. 1 we show the calculated photoluminescence intensity ratio /1//2 from the two wells as a function of the 
absorbed power density P for different values of the barrier width h. The self-consistent electric field F generated 
between the two wells in the same cases is shown in Fig. 2. When P is lower than a critical value, which depends on 
the barrier width, the photoluminescence intensity ratio /1//2 approximately decreases as P~^ and the electric field 
increases as P. When P reaches and exceeds the critical value an electric field of the order of 3 4 x 10^ V cm"-"^ gets 
established which slows down the tunneling of electrons from well 2 to well 1 and enhances the photoluminescence 
from well 2. At very high absorbed power density the ratio /1//2 tends to the value G1/G2 for two uncoupled wells. 
These results compare quite well with the experimental findings 1^ if one assumes a 0.1 % excitation efficiency. 

A deeper understanding of the behavior shown in Fig.s 1 and 2 can be reached by an approximate solution of the 
self-consistent procedure described above. Let us first concentrate on the dependence of the electric field, i.e. n, on 
the absorbed power density P. In the steady state we have ?T.2/''"e = P2/'''h which combined with Eq. (^) and Eq. 
(§d) gives 
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When P is not too high the second term in the squared root is small and we get the simple self-consistent equation 
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In a first approximation the smooth dependence of the tunneling rates on the electric field through the corresponding 
matrix elements can be neglected so that is constant while Tg still depends on n through the Fermi function. The 
argument of the Fermi function in Eq. (||) is negative until the alignment of the electron bands is achieved due to the 

electric field, i.e. in the region n < ricrit where e^b ricrit/^o^r = £'ei'' ~ E^j^ ■ In this range is approximately constant 
and much smaller than Th- On the other hand for n > Ucrit, the electron tunneling time increases exponentially with 
the adimensional parameter eFb/ksT, i.e. with n. The behavior of Te(n) imposes different solutions to Eq. (|lO|) at 
different values of the absorbed power density. At low P we have Te{n)P/hv ^ n and therefore n ~ ThP/hv. At high 
P we have Te{n)P/hv 3> n and therefore n ~ ricrit is constant. The critical value of the electric field is related to the 
critical value of the absorbed power density by Pcrit — hvricrit/Th- Since Th increases exponentially with the barrier 
width 6, Pcrit decreases exponentially with increasing 6 in agreement with Fig. 2. 

The behavior of n{P) allows us to understand also the features of Fig. 1. By approximating Eq. (^) with 

n2 ~ G2Te/ (1 -I- nX2Te) wc get 
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When P <^ Pcrit the electron tunneling time is almost constant, nX2Te <C 1 and the ratio /1//2 decreases as n~^, i.e. 
P~^. The behavior of /1//2 changes drastically at the critical power Pcrit due to the exponential change of Tc and, 
finally, for nA2Te 1 the asymptotic value /1//2 = G1/G2 is obtained. 
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For purpouse of comparison with the following Section, in Fig. 3 we show the behavior of the calculated normalized 
photoluminescence intensities /i//ioo, i=l,2, in the same ADQW of Fig.s 1 and 2 as a function of the barrier width b 
and for different absorbed powers. The normalization lioo is the photoluminescence current density ioi b ^ oo. The 
corresponding self-consistent electric field is shown in Fig. 4. The normalized photoluminescence intensity of well 2 - 
Fig. 3 - vanishes at smaller b when the tunneling current becomes higher. In this limit Je — > G2. At the same time 
the normalized photoluminescence intensity of well 1 tends to (Gi + Je) /Gi = I + G2/G1. The self-consistent electric 
field - Fig. 4 - gradually increases for decreasing b, until, in the tmmcling dominated limit /2//20C 0, it vanishes 
exponentially. 



III. TUNNELING FROM A QUANTUM WELL TO SURFACE STATES 

The considerations developed in the previous sections can be extended to describe the case where a single quantum 
well is close to a surface. This is another situation where tunneling followed by recombination takes place jf]-^,!). 
Electron and hole bands of defect states localized at the surface play the role of the el and hhl bands of the missing 
well. In a steady-state situation the system can be still described by the rate equations ([^) where the index 1 is 
associated with the surface and the index 2 with the well. Let ni and pi be the charge concentrations in the donor- 
like and acceptor-like surface bands. We can neglect the photogeneration of pairs at the surface so that Gi = 0. 
Moreover XiUiPi = Ii represents the nonradiative recombination current density at the surface. Since electrons and 
holes recombine (through a multiphonon process) very fast with respect to the other relevant time scales, one is 
allowed to take Xi —^ 00. 

The tunneling current densities from the occupied states of well 2 to the empty surface states 1 are still given by 
the general expressions (^) and (^). Now, however, the surface densities i^i^^ and are not constant in general. 
This implies the following definition for the corresponding Fermi energies: 
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for the donor-like band, and 
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for the acceptor-like band. As in the case of the two wells, we measure the energies from the bottom of the el or 
hhl band of well 2. More importantly, v'^'' and v\^^ depend on the material characteristics and can be orders of 

magnitude different. If <C v^l\ the effective mass difference between electrons and holes can be overcompensated 
by the difference in the surface densities and holes can tunnel more effectively than electrons do. This is the case, 
for instance, at the GaAs-oxide interface plp^]. In this case holes accumulate at the surface and in a steady-state 
situation pi > m. The electric field still given by Eq. (||) turns out to be negative. 

As in the case of the two wells, the tunneling current densities can be approximated by Je = ng/Te and = P2/Th- 
Now, however, we distinguish two cases. If ni > pi, the tunneling rates are given by Eq. (|) and Eq. (0) with 

ly^^^ = vi^^O) and — iylp{0). If pi > ni, we have 
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When the electric field and the tunneling rates are known, the solution of the initial rate equations is still given by 
Eq. (0). With the conditions Gi — Q and Ai ^ cx) we get the simpler formulas 
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P2 — n2 + n . (16d) 

The same recursive method explained in section II, allows us to find the steady-state values of the electric field, of 
the charge concentrations and of the tunneling rates. 

The evaluation of the luminescence intensity I2 implies the knowledge of the nature of the surface states. We can 
try to get information on the surface states by fitting experimental photoluminescence data. We will concentrate 
on the specific example of an Alo.sGaa.TAs surface with a nearby GaAs quantum well [QJ^. At energy close to the 
bottom of the el band of well 2 the Alo ^Gao iAs surface has only donor-like states belonging to the exponentially 
vanishing Urbach tail 



^ rni^^ ( AEc + eFb - EiV c, 
'^i'He) = ^exp{ ^ £i ) . (17) 



Such states are assumed to be nodal hydrogenic wavefunctions with radius Tg fixed by the depth into the gap. 
Their explicit expression is given in Appendix B. We assume that at the top of the gap the state density is the two 
dimensional density of free Alo.3GaQ,TAs electrons with effective mass mi^\ The parameter Ce will be considered as 
a fitting parameter. According to Eq. ([l7| ) and Eq. ([T^ ) where we can assume eI\^ — — cxd, the Fermi energy for the 
donor-like surface band is 

e^^^=AEe + eFb-Eif+eeln(^] . (18) 

\ml 'ee / 

On the other hand, at energy close to the bottom of the hhl of well 2 the AlQ.^GaojAs surface has a very high 
concentration of acceptor-like defect states [|l8|. We schematize them again by nodal hydrogenic wavefunctions |l9| ] 
but with radius to be considered as a second fitting parameter. These states are assumed to be distributed in 
energy with constant density vj^^ . The Fermi energy for the acceptor-like surface band is then 
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The tunneling matrix elements for the above surface states and surface densities are evaluated explicitly in Appendix 
B. 

According to the experiment reported in 0] we choose the well width a2 = 60 A the temperature T — 4.2 K, the 
photon energy hiy — 1.608 eV and the incident power density Pi = 0.5 Wcm^^. The incident cfhciency is estimated 
to be 1% and we take G2 = 0.01 Pi/hv. The relevant material parameters are Ai?c = 0.3 eV, AE^ — 0.128 eV, 
mi^^ = 0.091 TO, rrie^ = 0.067 to, ro'^'^ = 0.34 to, to being the free electron mass, and e,- = 12. Moreover we put 

A2 = 6 cm^ s" -1 0. We assume an acceptor-like surface state density = 10^'' cm" -2eV-i |2l| with E^^^^ ~ 1 eV 
(the results we found do not depend crucially on this particular value). The free parameters, and r^, are fixed by 
fitting the normalized photoluminescence intensity I2/I200 to the experimental data |Q] obtained for different values 
of the barrier width b. A least-square-error procedure gives the unique solution = 12 meV and — 11 A. In Fig. 
5, comparison is made between the ratio /2/-^2oo, calculated with these values, and the experimental data. 

Due to presence of two parameters the agreement between theoretical and experimental data in Fig. 5 should be 
considered as a source of information for these two parameters in the case the model is valid more than a proof of 
validity of the model itself. A reliable check of validity of our model is obtained by comparing the theoretical values 
of the self-consistently estimated electric field F with the values deduced from measured Stark shifts as a function of 
the QW excitation |^ . The calculated field values are shown in Fig. 6 as a function of the barrier thickness 6, for 
different values of the incident power. It is seen that, for high levels of excitation, the field approaches values of order 
10^ V cm~^ and keeps increasing when the barrier becomes thinner. Direct comparison with experiment is made for 
a 80-A-thick barrier, a choice dictated by the wide laser power range where data were available for this moderately 
tunneling sample. To obtain the value of the electric field from the measured Stark shift AEp we use the relationship 
empirically deduced from known results in the literature [p2| , p3| . The constant K is fixed by imposing 
that the electric field obtained from the measured Stark shift AEp — 0.33 meV, corresponding to an incident power 
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density Pi = 0.68 W cm ^, coincides with the calculated value. We obtain K = 1.25 • 10 ^ meV (V cm ^) which 
compares well with the values extracted from the data of Ref. if ~ 1 • 10~^ meV (V cm~^)~'^/^, and Ref. 
K -0.6- 10-'^ meV (V cm^i)-3/2. 

The comparison between the calculated and experimentally deduced electric field is shown in Fig. 7 as a function 
of the incident power density. The agreement is very good for a power range extending over 3 orders of magnitude. 
The near-coincidence of the scale factor K with other independently obtained values should not be attributed much 
importance, on one side because of the approximations used in the model, e.g. density and distribution of surface 
states, and, on the other side, because of the experimental Stark shift is observed in the luminescence from the 
quantum well, where the field is non- uniform and possibly different from its value in the barrier Instead, the 
correct functional dependence of the field on the excitation level is a clear confirmation that the self-consistent 
approach provides a reasonably accurate description of the whole process. 



The photoluminescence efficiency in asymmetric double quantum wells and in near-surface quantum wells is strongly 
influenced by the tunneling of both electrons and holes between the two wells, or the well and the surface states. 
The theoretical model presented here takes into account this effect and allows a quantitative prediction for the 
photoluminescence rates and the value of the electric field which needs to be established across the tunneling barrier 
in a steady-state situation. Experimental results regarding both ADQW and near-surface QWs, namely changes in 
photoluminescence efficiency and peak shifts due to the self-induced electric field, are very well accounted for. 

Our model can be readily generalized to include the case of an externally applied electric field as well as the case 
of not constant excitation. 



We call z the coordinate orthogonal to the interfaces and suppose that the quantum well 1 is in —ai < z < and 
the quantum well 2 in b < z < b + a2- Firstly let us consider the case of electrons. As schematized in Fig. 8 we 
assume a rectangular potential profile with left and right discontinuities VJ*-^-* = AEc and Vr^'' = Ai?c — eFb/2 for 

(2) (2) 

the well 1 and = AEc + eFb/2 and Vr — AEc for the well 2 where F is the electric field in the barrier region 
< z < b. The electron wavefunctions at energy e, measured from the bottom of the el band of well 2, are 



IV. CONCLUSIONS 
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APPENDIX A: ADQW MATRIX ELEMENTS 
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for the well 1 and 
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for the well 2. Energy iJ^*^'' measures the bottom of the el band in the well z = 1, 2 from the bottom of the same well 
and is determined by solving 
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The phase shifts are J'-*-' — tan ^{k^"^^ /kf^), i ~ 1,2. The constants C'-'\ « = 1,2, are fixed by normalizing the 



wavefunctions 
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For e = and with the assumption that in Eq. (j^) vanishes everywhere but in the barrier region where Vg = Ai?c 
the tunneling matrix element is 

(<i>«|T4|$[,'^) = 
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Jo Jo VA VA 



A 



3in[fc(i)ai +(5(i']sinj(2) 



^(2) „ ^(1) 



ifc. 



(1) 



(A13) 



The last two factors in this expression quickly oscillates with fci^"* 
,(1)^1.(1)^ _ J2rne\Ef^ - £;^P - eFb)/n^ - (fci^V- By taking the value for ki^^ = (maximum value) the 
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electron tunneling matrix element is estimated as 
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In the case of holes we have a completely analogous situation where the relevant band is hhl instead of el. The 



expressions derived for electrons still hold for holes with the substitutions m 



eFb 



^eFb. 



1,2, AEc 



AK, 
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APPENDIX B: NEAR-SURFACE QW MATRIX ELEMENTS 



Following the notations of Appendix A, the surface is defined by the plane z — and the well 2 is in 6 < z < 6 + a2. 
Firstly, we consider the case of electrons and we measure the energy from the bottom of the el band of well 2. The 
state ^i^\x,y,z) is given by Eq. ( [A6| ). The donor-like surface state ^i^\x,y,z) is approximated by a truncated 2p 
hydrogenic wavefunction 



5/2 



z <0 

1 0< z 



(Bl) 



where r ~ \/1e^^^^^]P^^V^ ■ The state is at energy h? /{^rrie rD below the bottom of the conduction band for the 

barrier material where the electron effective mass is rn'e^ . By imposing this energy to correspond to e we determine 
the radius 



Srrie' ( ^Ec - E^t' + eFb - e 



(2) 



(B2) 



Using parabolic coordinates r] = (r — z)/re, ^ ~ (r + z)/re, = tan ^{y/x), the tunneling matrix element between 
the well and surface states at e = is: 



r)+2b/re 

drj 



d£_ / dip 



rl{^ + rj) (^_^)e-(?+'?)/4 ^(2) g^^^ j(2)gfeP'[«-^K/2-6] 

iaE, 



Q r- 3/2 



r^AEcC^^') sin (5(2) 



e ' 



1 



26/re 



4(2fcf Ve - 1) 
b/re 



32(2/cf Ve - If 
1 



1 



4(2fc, 



1)2 32(2fcpVe - 1)3 



(B3) 



In the case of holes we have a completely analogous situation where the relevant band of well 2 is hhl instead of 
el and the acceptor-like surface state is given by Eq. (Bl) with —^ r^. Equation (B3) gives the tunneling matrix 
element for holes with the substitutions re ^ rt, AEc AEy, eFb —eFb. 
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FIG. 1. Calculated emission intensity ratio I1/I2 from two coupled asymmetric wells vs absorbed power density P for various 
values of the barrier width b. 



FIG. 2. Calculated electric field F which establishes between the two wells, in the same cases of Fig. 1. 

FIG. 3. Calculated normalized photoluminescence intensities /1//100 and I2/I200 vs barrier thickness b for different absorbed 
power densities P. Uoo is the photoluminescence intensity from well i in the limit 6 ^ 00. 



FIG. 4. Calculated electric field F which establishes between the two wells in the same cases of Fig. 3. 



FIG. 5. Normalized photoluminescence ratio I2/I200 of a near-surface well vs the surface-barrier thickness b. Dots: experi- 
mental data from Ref. [1]; solid line: best fitting in terms of the self-consistent model. Incident power density is Pi = 0.5 W 
cm~^. 

FIG. 6. Calculated electric field F across the surface barrier vs the surface barrier thickness b for different incident power 
densities Pi. 



FIG. 7. Comparison between the electric field F calculated from the model (solid line) and deduced from the measured Stark 
shifts [8] (dots) vs the incident power density Pi. The surface-barrier thickness is b = 80A. 



FIG. 8. Energy-space diagram for electrons in the asymmetric double quantum well. 
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